NON-MARKOV PROPERTY OF CERTAIN EIGENVALUE 
PROCESSES ANALOGOUS TO DYSON'S MODEL 
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Abstract. It is proven that the eigenvalue process of Dyson's random matrix 
process of size two becomes non-Markov if the common coefficient l/v2 in the 
non-diagonal entries is replaced by a different positive number. 
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1. Introduction 

Dyson [3] has introduced the matrix-valued stochastic process 
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to model the dynamics of particles with the Coulomb type interactions, where Bi/s 
are real Brownian motions and B^j's for i < j are complex Brownian motions all of 
which are mutually independent. He proved that the eigenvalue processes Aj, . . . , Ajv 
satisfy the (system of) stochastic differential equations 



dm + 2 



2 ^ Ai(t) - \j(t) 



dt 



with (3 = 2. It has been proven later that if the complex Brownian motions are 
replaced by real or quaternion Brownian motions, the eigenvalue processes satisfy 
similar stochastic differential equations with (3 = 1 or 4, respectively. (See [H H] 
for discussions based on the stochastic analysis.) These processes are now called 
Dyson's Brownian motion models for GOE, GUE, and GSE when j3 = 1,2, and 
4, respectively. In any case, it is remarkable that the process A — (Ai, • • • , Ajy) is 
Markov. 

We may ask the following question: "Does the process A remain Markov if we 
replace the common coefficient l/y/2 by a different positive number?" In this paper, 
we give the negative answer to this question when the matrix size N = 2. 

Let c > and 8 > 0. Consider the 2 x 2-matrix-valued process 
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where B x and B 2 are two independent standard Brownian motions and £ s is a Bessel 
process of dimension 6 starting from which is independent of B x and B 2 . We see 
in Lemma [2.21 that S c ' 5 with 6 = 1, 2, or 4 is unitarily equivalent in law to 



c , Su ,_l B x {t) Jc/2B 3 (t) 



;i ' 2) HC ' (t) 1 'c/2fl s (f) 5aW 

with i?3 a real, complex, or quaternion Brownian motion independent of B\ and _B 2 , 
respectively. Let \i(t) and \2{t) for t > denote the eigenvalues of the Hermitian 
matrix E c,s (t) such that Xi(t) > A 2 (t). Define the two-dimensional process A c,<5 = 
(Ai,A 2 ). 

When c = 0, X\(t) and A 2 (t) are nothing but the order statistics of B x (t) and 
B 2 {t), that is, Ai(t) = max {B x (t), B 2 (t)} and A 2 (t) = min {B x (t), B 2 (t)}. Hence it 
is obvious that the process A ' 5 is Markov. 

When c = 1, the process (11.11) is a time- dependent version of Dumitriu-Edelman's 
matrix model for beta-ensembles (cf. [2]) and we see in Lemma [2TT1 that the processes 
Ai(i) and A 2 (t) satisfy Dyson's stochastic differential equations with index (3 = 6 
given by 

(1.3) dX ^= d ^ + 2{Ut)-Ut)) dt ' 

(1.4) rfA2W ^ W + _j_ df 

for two independent Brownian motions (3\{t) and /3 2 (t). In particular, the process 
A li5 (t) is Markov. 

Theorem 1.1. The process A c ' 5 is Markov if and only if c G {0, 1}. 

We prove this theorem by reducing it to the following. 

Theorem 1.2. Let 61,62 > 0. Let X Sl and Y &2 be two independent squared Bessel 
processes starting from of dimension 6\ and 62, respectively. Then the process 
Z c (t) = cX Sl (t) + Y S2 (t) forc>0 is Markov if and only if c G {0, 1}. 

Theorems 11.11 and 11.21 seem similar to Matsumoto-Ogura's cM — X theorem [6]. 
Let X be a Brownian motion and set M{t) = sup 0<s<t X(s). When c = 0,1,2, the 
process cM — X is Markov; indeed, —X is a Brownian motion, M — X is a reflecting 
Brownian motion by Levy's theorem (see, e.g., [7J Thm.VI.2.3]), and 2M — X is a 
three-dimensional Bessel process by Pitman's theorem (see, e.g., [7J Thm.VI.3.5]). 

Theorem 1.3 (|6J). The process cM — X is Markov if and only if c G {0, 1, 2}. 

2. Non-Markov property of the eigenvalue processes 

Proof of Theorem li.il provided Theorem li.^l is justified. An elementary calculation 
shows that Ai and A 2 are given by 

Ai(t) =\ {Bi(t) + B 2 (t) + y /(B 1 (t)-B2(t)y + 2c^(ty} , 
A 2 (t) =- {B 1 (t) + B 2 (t) - v /(B 1 (t)-B 2 (t)y + 2^(ty} . 
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Set B 3 (t) = {B 1 (t) + B 2 (t)}/V2, X^t) = {B 1 (t) - B 2 (t)} 2 /2 and Y*(t) = ^(t) 2 . 
Then B% is a real Brownian motion, X 1 and Y s are squared Bessel processes of di- 
mension 1 and 5, respectively. Moreover, B%, X 1 , and Y s are mutually independent. 
It follows that 



X 1 (t) =±= |5 3 + y/x^+cY'it)} 
X 2 (t) =-L |5 3 _ y/ X i(t) + cY'(t)} 



It is obvious that the two dimensional process A c > s = (Ai,A 2 ) is Markov if and 
only if so is the process (Ai + A2, Ai — A2). Since 

(2.1) X 1 + X 2 =V2B 3 , 

(2.2) Ai - A 2 =v / 2 VX 1 + cY s 

and they are independent, for the process A c ' s to be Markov it is necessary and 
sufficient that the process X 1 + cY s is Markov. This is equivalent to c = or 1 by 
Theorem 11.21 □ 

Lemma 2.1. For c = 1 and 5 > 0, consider the 2 x 2 -matrix-valued process H 1 '* 
defined by (11. ip . TTien i/ie corresponding eigenvalue processes satisfy the stochastic 
differential equations (II. 3p - (11.41) . 

Proof. Set A = (Ai — X 2 )/V2. Then, by (12. 2p for c = 1 and by Shiga- Watanabe's 
theorem (see, e.g., [71 Thm.XI.1.2]), we see that the process A is a Bessel process of 
dimension 1 + 5. Hence we have 

(2.3) dX(t) = dB 4 (t) + -J—dt 

2 X(t) 

where B± is a real Brownian motion independent of B3. If we set f3% = (B3 + B4) /\/2 
and (3 2 = (B 3 — B 4 )/y2, then (3\ and f3 2 are two independent real Brownian motions. 
Therefore, combining (12.31) with (12. ip . we conclude that (ll.3p - (ll.4l) hold. □ 

Lemma 2.2. Let c > 0, 5 = 1, 2, or 4, and S c ' 5 and S c,<5 be the matrix-valued pro- 
cesses defined by (11. ip and ( 11. 2ft . respectively. Then, there exists a unitary matrix- 
valued process U$(t) such that 

(~^(t)) ^ (u s m^(t)u* 8 (t)) . 

V / t>o V / t>o 

In particular, eigenvalue processes associated with E C ' S and E c,s have the same law. 
Proof. We define 

Us(t) = ^ B 3 (t) J 1b 3 (*)^o + (^q 1 J ls 3 (*)=o 



[Bs(t)| 

by using .63 in (II. 2p . Then we have 



B 1 {t) ^fc~/2\B z (t)\ 
G/2\B 3 {t)\ B 2 (t) 



U 8 m^(t)U* s (t) = 

which shows the desired result since \B 3 \ L = £ d . □ 
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3. Transition probability density of squared Bessel processes 

In this section, we recall some basic asymptotic estimates on the transition prob- 
ability density pf(x, y) of squared Bessel processes of dimension 8 which we shall use 
later. We first note that it has an expression 

1 /yx^- 2 )/ 4 / x + y\ ( ' y/xy 



for x, y > 0, where l v stands for the modified Bessel function of index v (see, e.g., [H 
Cor. XI. 1.4]). Now let us recall following two asymptotic estimates on the modified 
Bessel function (see, e.g., Sect. 5.16.4 of [5]): 

(3-2) Iv{x) ~ T(v+1) (^Y aSXi °' 

(3.3) I u {x) ~ —j= as x t oo. 

\j2~KX 

Here, f(x) ~ g(x) means f(x)/g(x) — > 1 in the subsequently indicated limit. 
Using (13.21) in (13. ip . we can derive 

< X4 > rf(° + -»> = (a )W(f/2) B[p (-5 

for t,y > and 
(3.5) 



lim y l - s/2 p & t (x,y) = x l ~ s/2 p s t (0+, x) 



a; \ 

exp l I 



(2t) tf /2r(5/2) ^ V 2tJ 
for t,x > 0. On the other hand fl3.3H together with ( 13. ip yields 

,, R s n 1 V (S ~ m ( x + y-2^xy\ 



as y'xy — >■ oo. 

4. Non-Markov property of weighted sums of two independent 

squared Bessel processes 



For the proof of Theorem 11.21 we may restrict ourselves to < c < 1; other- 
wise consider Z c /c instead. We prove that Z c is non-Markov by checking that the 
conditional law 

(4.1) P (Z c (2) G dz 3 | Z c (e) = z x , Z c (l) = z 2 ) for < e < 1 

does depend on (e, z\). This conditional law has the density 

P{Z\2) G dz 3 | Z\e) = Zl , Z\l) = z 2 ) = q{ J> Z3;£ >* l] dz 3 , 

where q(z 2 , z 3 ; e, Z\) and q{z 2 \ e, z-y) are the densities of the joint laws of (Z c (s), Z c (l), Z°(2)) 

and (Z c (e), Z c (l)), respectively. Thus it suffices to prove that the fraction q(z 2 , z 3 ; e, zi)/q(z 2 ; e, z\) 

depends on (e, z\). 
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To this end, we shall use the integral expression 

rz\ pz2 pzs 

q(z 2 ,z 3 ;e,z 1 ) = / dx x \ dx 2 / dx 3 A 1A A li2 A li3 , 



q(z 2 ;s,zi) = / dxi / dx 2 A ltl A lt2 
Jo Jo 



where 



A 1A = p ^(0 + , Xl )p^(0 + , Zl - cxx), 

A lt2 =p[^ s (xi,x 2 )p^_ £ (zi - cxi, z 2 - cx 2 ), 

A lt3 =p 5 l 1 (x2,x 3 )p 5 l 2 (z 2 - cx 2 , z 3 - cx 3 ). 

We divide the proof into several steps. First of all, we prove 

Lemma 4.1. Let /(A,-) for A > be a bounded measurable function on (0,1). 
Suppose that f (A, x / A) converges to a constant f (00, 0) for any x G (0, 1) as X — > 00. 
Let <p G C l ((0, 1)) and suppose that 0(0+) = a Gl, <p'{0+) = b > and <p'(x) > 
for x G (0, 1). Let v > 0. Then 

e- x ^f(X, x)x v - l dx ~ /(oo, 0)^A~ ly e- aA as A -> 00. 

Proof. Changing variables to u = Ax, we find that the left hand side of (14. 2 p equals 

X~ u e- aX [ e- XWu/x) - a} f(X,u/X)du. 



Note that A{0(w/A) —a}> Ku for u G (0, A) and A > where K = inf xe ( Oil ){0(x) — 
0(0+) }/x > 0. Hence we see that 

/* A /»oo 

lim / e~ x{ ^ u/x) ~ a] f{X,u/X)du = /(oo,0) / e^u" -1 ^ 

by the dominated convergence theorem. □ 

Second, we take the limit as e — > 0. 
Lemma 4.2. 

g(32,z 3 ;e,2i) 9(22, 23; *i) 

hm — — = — r— 

£^0+ q(z 2 ;e,z 1 ) q{z 2 ;z 1 ) 

with 

q(z 2 ,z 3 ]z 1 )= I dx 2 I dx 3 A 2A A 2y2 , q(z 2 \z x ) = I dx 2 A 2>1 
Jo Jo Jo 

where A 2t2 = A± j3 and 



A 2 ,i = A lt2 
Proof. We know that 



= Pi 1 (0+,X2)pi 2 (zi, z 2 - cx 2 ). 

e^0+,x^0+ 



(x 1 )^- 1 (z 1 ~cx 1 Y s ^- 1 ( 1 r n , . 

Al > 1 = (^(wrMrw exp l"2i {Z1 + (1 ~ c)xi} 
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from (13 .4p . Now we can rewrite q(z 2 , z 3 ; e, Zi)/q(z 2 ; e, Zi) as Fi/G\ with 

(4.3) F x = I ' A M (e, x l )x ( f l/2yi e^ /£)xi dx l 

Jo 

(4.4) G 1 = r A 1>5 (e, x 1 )x[ Sl/2) ~ 1 e- { " /e)xi dx 1 

Jo 

where c= (1 — c)/2 and 

Ai i4 (e,xi) =(21 - cxi) (<52/2)_1 / dx 2 dx 3 Ai >2 A lj3 , 

Jo Jo 

1>Z2 

A 1>5 (e,xi) =(zi - cxO^- 1 / <*Mi,2- 



Using Lemma 14.11 in the integrals H4.3f) and (14.41) , we have 

F 1 ~e* l/2 r(<fi/2)r* /2 A M (0, 0), 

d ~ £ < 5 i/ 2 r(5 1 /2)c- 5l / 2 A li5 (0, 0) 

as e — > 0+. Here we have used the fact that Ai^(e, Xi) and A 15 (e, xi) are continuous 
in e G [0, oo) and x\ G [0, zi]. Therefore, F\/G\ approaches to A 1)4 (0, 0)/A lj5 (0, 0) = 
q{z2 i z 3 ;zi)/q(z2;z 1 ). □ 

Third, we study the asymptotic behavior of the numerator q(z 2 , z 3 ; z\) as z 3 — ► 0+. 

Lemma 4.3. 

lim Z3~ (<5l+<52)/2 g(2;2, 2: 3 ; zi) = Cife; Zi) 
Ci = I >u^^-\l - cufz^du, q{z 2 ; z x ) = [* dx 2 A 3 , 1 A 3>2 



o 



where A 3 ^\ = A 2t i and 

A 3 , 2 = {x 2 ) 1 - 5 ^\z 2 - cx 2 ) 1 -^ 2 p 6 1 ^0+,x 2 )p s 1 *(0+,z 2 - cx 2 ). 
Proof. Recall that 

(4.5) q(z 2 , z 3 ; z x ) = dx 3 A 2t3 (z 3 , x 3 ) 

Jo 

where 

^2,3(^3,^3)= / dx 2 A 3A p s 1 1 (x 2 ,x 3 )p s 1 2 (z 2 - cx 2 ,z 3 - cx 3 ). 
Jo 

Here we note that ^1 does not depend on z 3 nor x 3 . If we take x 3 = z 3 u for 
< u < 1, we have 



A 2 , 3 (z 3 , z 3 u) = / dx 2 A 3A p s 1 1 (x 2 ,z 3 u)p 5 1 2 (z 2 - cx 2 ,z 3 (l - cu)). 
Jo 

Using (13.51) . we have, as z 3 — > 0+, 

,3 2 -^)/% 3 (,3, Z 3 U) - ^^(l - CU)^- 1 / ^^3,^3,2. 
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Changing variables to u = x 3 /z 3 in the integral (14.51) . we obtain 

Jo 

which converges to Ciq(z 2 ; Z\) as z 3 — > 0+. □ 

Fourth, we study the asymptotic behaviors of q{z 2 \ z\) and q{z 2 \ z\) as z 2 — > oo. 
Recall that 

q(z 2 ; zi) = / dx 2 A 3jl A 3t2 
Jo 

= r dx 2 x l 2 - Sll2 {z 2 - cxs) 1 -^! 1 (0+, x 2 ) 
Jo 

x p s 1 2 (z 1 ,z 2 - cx 2 )p 5 1 1 (0+,x 2 )p s 1 2 (0+,z 2 - cx 2 ) 

3-(5i+5 2 )/2 f 1 r l-Si/2 /i \l-8 2 /2 Sun, \ 

=z 2 "I duu l/ (1 — cu) 1 Pi(0+, z 2 u) 
Jo 

x p{ 2 (z h z 2 {l - cu))p{ 1 (0+, z 2 u)p{ 2 (0+, z 2 {l - cu)) 

and that 

q(z2] zi) =z 2 / dup{ 1 (0+, z 2 u)p{ 2 (zi, z 2 (l - cu)). 
Jo 

Lemma 4.4. Let r > 0. Then 

1 A a\ Q\ z 2'i Z 2 r) \-5i/2 -z 2 /2 

(4.6) — ~ C 2 D{r) 11 e 2/ as z 2 00 

g(z 2 ; z 2 r) 

where C 2 is some positive constant depending only on Si and 62 and 

1 - c 



D(r) 



1 — c + ■sfrc 
Proof. If we express q(z 2 ; z 2 r) as 

r (l- 52 )/4^-l)/2 f 1 fl ^ u)e -^Mu) u Si/2-l du 
JO 

using 

0i (w) = &i« + \/r {l — y/l — cm} + ai 

with b\ = \ — c and a! = (-^/r — l) 2 /2 + 1/2, then fi(z 2 , ■ ) turns out to be a bounded 
continuous function such that fi(z 2 ,u/z 2 ) converges to a constant depending only 
on Si and S 2 as z 2 — > 00, by (13. 6p . Since 0i and /1 satisfies the assumptions, we can 
use Lemma 14.11 and hence we obtain 

(4.7) q{z 2 - z 2 r) ~ C 2 ,ir il - 52)/4 <j)' 1 (Q+)- 5l/2 z 2 1/2 e- aiZ2 as 2 2 -> 00 

with some constant C 2) i depending only on Si and S 2 . 
We also have a similar expression 

r (l-* 2 )/4^-l)/2 u)e - Z2Mu)u S l/2 -l du 
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for q(z 2 ; z 2 r) using 



4> 2 (u) = b 2 u + \/r {l - \Jl- cu] + a 2 

with b 2 — (1 — c)/2 and a 2 = (y/r — l) 2 /2 and a function 72(^2, • ) as before. Thus 
the same argument yields 

(4.8) q{z 2 - z 2 r) ~ C 2 ,2r (1 - 52)/ y 2 (0+)- 5l/2 z~ V V a2 * 2 as 2 2 -> 00 

with some constant C 2)2 depending only on 5\ and 5 2 . 

Using (g2D and (jM])' together with 0i(O+) = 61 + ^/rc/2 and 0' 2 (O+) = 6 2 + v/rc/2, 
we obtain (|£Bj) . □ 

Now we are in a position to prove Theorem 11.21 

Proof of Theorem ! 1.21 . Let < c < 1. We combine Lemmas I4.2[ 14.31 and 14.41 to 
obtain 

lim e^ 2 lim ^+*)/ a lim g^ifl^j = C 3 L>(r)-^ 2 
z 2 ^oo 23-^0+ e-*o+ g(^ 2 ; e, z 2 r) 

for some constant C3 which depends only on 5\, 5 2 and c. Therefore we conclude 
that the conditional probability (14. ip does depend on (e, z\), which proves that Z c 
is non-Markov. □ 
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